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ABSTRACT 

We investigate the dynamical effects of dark matter subhalos on the structure and evolution of 
a galactic disk, using semi-analytic method that includes approximated and empirical relations as 
achieved in detailed numerical simulations of the cold dark matter model. We calculate the upper 
limit for the size of a galactic disk at a specific redshift z, based on the orbital properties of subhalos 
characterized by their pericentric distances from the center of a host halo. We find that this possibly 
largest size of a disk as determined by the smallest pericentric distances of subhalos shows the charac- 
teristic properties, which are basically in agreement with an observed galactic disk at low and high z. 
Namely, it is found that a massive disk can have a larger size than a less massive one, because of its 
stability against the destruction effect of subhalos. Also, with fixed mass, the size of a galactic disk at 
low z can be larger than that at high z, refiecting the orbital evolution of subhalos with respect to a 
host halo. These results suggest that the presence and structure of a galactic disk may be dynamically 
limited by the interaction with dark matter substructures, especially at high z. 

Subject headings: cosmology: dark matter — galaxies: formation — galaxies: structure — galaxies: 
interactions 



1. INTRODUCTION 

The cold dark matter (CDM) model is now a standard 
theory for the structure formation of the universe. The 
model provides us with a basic theoretical framework for 
understanding the hierarchical clustering of dark matter, 
where large dark halos are assembled via merging and 
accretion of smaller halos. Indeed, the prediction of the 
CDM model is in excellent agreement with observations 
of large-scale structures (Tegmark et al. 2004). How- 
ever, recent high-resolution N-body simulations based 
on the CDM model have highlighted that as a result of 
these merging processes, a large dark halo like that of 
the Milky Way contains numerous dark matter substruc- 
tures or subhalos inside its virial radius, in contrast to 
the small number of known Milky Way satellites (Klypin 
et al. 1999; Moore et al. 1999). This so-called missing 
satellite problem remains unsolved yet, although various 
ideas are proposed to solve it based on, e.g., astrophysi- 
cal baryonic processes such as UV feedback to suppress 
star formation and/or observational selection effects for 
faint satellite galaxies in the Milky Way (e.g.. Bullock, 
Kravtsov & Weinberg 2000; Madau et al. 2008; ToUerud 
et al. 2008; Koposov et al. 2008; Koposov et al. 2009; 
Maccio et al. 2009). 

If such dark matter subhalos indeed exist around a 
disk galaxy and some of these, especially massive ones, 
interact with a galactic disk in the course of their or- 
bital motions, then the stellar component of a galactic 
disk can be made so thick due to dynamical heating that 
it no longer exists as a thin stellar disk like that of the 
Milky Way at the present epoch. This issue has been in- 
vestigated by many researchers as a constraint on models 
of galaxy formation and evolution (e.g., Toth & Ostriker 
1992; Font et al. 2001; Ardi et al. 2003; Hayashi & Chiba 
2006; Villalobos & Helmi 2008; Hopkins et al. 2008; Pur- 
cell et al. 2009). Hopkins et al. (2008) recently showed, 
based on the realistic (radial) orbits of subhalos as ob- 



tained from cosmological simulations, that disk heating 
is less powerful than previously thought by Toth & Os- 
triker (1992), which is based on the simple assumption 
of rigid satellites with circular orbits. To test what this 
result implies for satellite accretion with high mass. Pur- 
cell et al. (2009) made high-resolution numerical sim- 
ulations for satellite-disk interactions and showed that 
in reality accretion events of mass ratio ^1:10 do not 
preserve thin disk components. Indeed, Stewart et al. 
(2008) demonstrated, also based on N-body simulations, 
that a majority 70%) of galaxy-sized halos with mass 
M ~ IO^^Mq at z = have accreted at least one ob- 
ject with M > IO^Mq over the last 10 Gyr and 95% 
have accreted an object with mass greater than Milky 
Way disk (M > 5 x IO^^Mq). As these works suggest, 
a galactic disk is a fragile system, so that its dynami- 
cal evolution may be severely limited by these subhalos. 
In particular, the presence of a thin galactic disk at the 
present epoch implies the absence of the effect of subha- 
los within its size, which suggests in turn that the size 
of a galactic disk may be regulated by the orbital radii 
of subhalos. Therefore it is important to investigate the 
relation between the size evolution of a galactic disk and 
dynamics of subhalos within a host halo. 

Several models for the structural evolution of a galactic 
disk have been studied (Mo, Mao & White 1998, here- 
after MMW98; Bouwens & Silk 2002, hereafter BS02; 
Kampakoglou & Silk 2007). The evolution models can be 
divided into two alternative approaches, one is forward 
approach and the other is backward approach. The for- 
ward approach is based on the CDM model, where both 
collapsing dark matter and baryonic gas acquire the an- 
gular momentum through tidal force and mergers. When 
a gas component cools, condenses and forms stars, then 
the system's angular momentum halts the collapse, lead- 
ing to the formation of a rotationally supported disk (Fall 
& Efstathiou 1980). Under the assumption that the frac- 
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tions of disk mass and angular momentum in the disk rel- 
ative to the halo, together with the spin parameter of the 
halo, do not vary during this collapsing phase. Mo, Mao 
& White (1999) obtained the following relations for the 
scale length of an exponential stellar disk: Rd oc H{z)~^ 
at a fixed circular velocity or Rd oc H{z)~^^'^ at a fixed 
halo mass, where H(z) is a Hubble parameter at z. On 
the other hand, the backward approach is based on the 
detailed models of the properties of local disk galaxies, 
e.g., gaseous, stellar, and mctallicity profiles, as well as 
current star formation rate (SFR) and age-metallicity re- 
lationship. This approach uses the local universe as a 
reference. Then the properties of a galactic disk at high 
z are derived on the basis of those at present time and 
backward calculation in time. Using the Milky Way as 
reference and assuming the continuous infall of metal-free 
gas from outside, BS02 provide the following size-redshift 
relation for a disk: r(z)/r(0) = 1 — 0.27 z. 

Observational studies of a galactic disk have also been 
put forward by several researchers. Using the Sloan Dig- 
ital Sky Survey (SDSS), Shen et al. (2003) showed the 
local size-mass relation, where a galactic disk with large 
mass has a larger scale length than that with small mass. 
Based on the Galaxy Evolution from Morphology and 
SEDs (GEMS) survey, Barden et al. (2005) presented 
the size-mass relation up to z ~ 1. Also, using the Faint 
Infrared Extragalactic Survey (FIRES), Trujillo et al. 
(2004) measured the size-mass relation in the rest-frame 
optical up to z ~ 2.5. Trujillo et al. (2006) presented 
that at a given mass the mean size of a galactic disk was 
~ 2 times smaller at z 2.5 than that we see today, 
using the result of FIRES, GEMS, and SDSS. 

It is worth noting that models in both forward and 
backward approaches are in good agreement with these 
observational results. However, there are a couple of un- 
resolved issues in these models. In forward approach, the 
model assumes that a gas disk forms a stellar disk instan- 
taneously at the epoch of a halo virialization, whereby 
only the evolution of a host halo is taken into account. 
In backward approach, the model docs not include the 
merging history of a host halo. Thus in these models, the 
important dynamical effects of subhalos on the structural 
evolution of a galactic disk are unclear. In particular, 
since subhalos are able to destroy a galactic disk and 
the probability of such events is expected to be high, the 
previous models are yet insufficient to understand the 
evolution of a galactic disk. Therefore, it is of great im- 
portance to fully take into account the effect of subhalos 
on the evolution model of a galactic disk. 

We remark here that this aspect of a disk evolution is 
automatically included in high-resolution gasdynamical 
simulations of a collapsing galaxy in the context of the 
CDM model. However, many of such simulation works 
have produced a disk which is too small compared with 
a galactic disk at present time (e.g., Navarro, Frenk & 
White 1995), most probably because the simulation re- 
sults are largely affected by yet poorly understood physi- 
cal processes, such as star formation and supernova feed- 
back, that work below the resolution limit. Indeed, Gov- 
ernato et al. (2007) showed that their simulation taking 
into account improved feedback models along with high 
numerical resolution largely minimizes the issue of a too 
small simulated disk. Here, to avoid such complexity 



in understanding the calculated results, we concentrate 
on the dynamical effect of subhalos on a galactic disk 
and adopt a semi-analytic method that includes approx- 
imated and empirical relations as achieved in detailed 
numerical simulations of the GDM model. 

The outline of this paper is as follows. In ^we describe 
our semi-analytical model that we use for the merging 
and accretion of dark matter halos. In ^we present the 
examples of our semi-analytical model and the results 
for the disk evolution. We summarize and discuss the 
implications of the current results in 21 

2. MODEL 

We describe here the properties of subhalos in a host 
dark halo, in particular the mass accretion history and 
the orbital evolution of subhalos after they enter the 
virial radius of a host halo. For our current work, 
we adopt the semi-analytic model developed by Zent- 
ncr & Bullock (2003), hereafter ZB03, with some mod- 
ification as explained below. In the foUowings, we 
adopt the standard set of cosmological parameters: the 
density parameter at present f2mO = 0.3, cosmological 
constant at present A = 0.7, Hubble constant Ho = 
lOO/i km s^^Mpc^^ with h = 0.72, baryon fraction at 
present fl^h^ = 0.02, and rms density variation averaged 
over 8h~^ Mpc as as = 0.95. 

2.1. The ZB03 Model 

2.1.1. Construction of Merger Trees 

Following ZB03, we track the mass accretion history 
of a host halo using the extended Press-Schechtcr (EPS) 
theory (Bond et al. 1991; Lacey & Cole 1993, here- 
after LG93) and adopt the merger tree algorithm pro- 
posed by Somerville & KoUat (1999, hereafter SK99). 
Using this method, we generate a list of the masses and 
accretion redshifts of all subhalos that have merged to 
form a host halo of a given mass at a given redshift. In 
the EPS formalism, the linear density field 5{M) as a 
function of the smoothing mass M is regarded as a ran- 
dom Gaussian field, which is specified by the mass vari- 
ance S{M). To derive this, we consider a density field 
smoothed with a spherical top-hat function with a radius 
R. Also, according to the spherical collapse model, the 
density field collapses and forms a virialized object at z 
once the density field 6 exceeds lo ~ (5crit(^) — 1-68 (see 
LC93). Then, denoting AS" S{M) - S{M + AM) and 
Slu = uj{t) — uj{t + At), the probability that a halo having 
mass M at time t has accreted the mass corresponding 
to a step of AS" in a time step associated with Su; is given 
as 



PiAS,SLu)d{AS) 



6lu 



27rAS'3/2 



exp 



(Slo) 



2AS 



diAS). 



(1) 

Following SK99, we generate merger trees by starting 
the redshift 2; = at which we consider the final halo 
and step backward in time. To reproduce the prediction 
of the conditional mass function of the EPS model, we 
must choose the appropriate time step. In this model we 
choose the time step given by Taylor & Babul (2004), 
Slu = [0.21ogjQ (M/Minin) + 0.1]Sljo, where M^in is a 
resolution limit mass and o^o = {dS / dM)\M Maihi- At 
each timestep we take a mass S{Mp) ~ S{M) + AS from 
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the probability equation (P). We treat events with A/p < 
Mmin as accreted mass and retain all information about 
mergers with Mp > Mmin- In this way, we generate the 
list of progenitor masses and accretion redshifts at each 
time step. This process continues until the masses of all 
progenitors arc less than Mmin- At each time step wc 
identify the most massive progenitor with the host halo 
and all less massive progenitor with accreted subhalos. 

We note here that the recent work by Zhang, Fakhouri 
& Ma (2008) showed cautionary remarks in the use of 
the SK99 model, which overestimates the abundances of 
small progenitor halos as large as a factor of about 2 com- 
pared to the prediction of the EPS formalism. Although 
a quantitative estimation for the effect of adopting dif- 
ferent Monte Carlo algorithms from SK99 is beyond the 
scope of this work, our results based on the distribution 
of pericentric orbital radii of subhalos, as detailed later, 
would not be sensitive to their precise abundances at 
each rcdshift. 

2.1.2. Density Distribution of Halos 

The spherical collapse model provides the density of a 
virializcd region for a dark halo. The density of a viri- 
alized halo is given as pvir = Pm^viv, where pm is the 
mean matter density of the universe. For Avir, we use 
the fitting function provided by Bryan & Norman (1998), 
where for our cosmological parameters, Avir ~ 370 at 
z = and Avir ~ 179 at high z. Then the virial radius 
of a virializcd halo with mass Mvir at z is 



i?v 



(3A/vir/47rAvirPm)^/^ 



(2) 



and the circular velocity at the virial radius, so-called 
virial velocity, is Kir = (GA/vir/i?vir)^^^- 

Recent numerical simulations discover the several an- 
alytical density profile of a dark halo (Navarro, Frenk 
& White 1997, hereafter NEW; Moore et al. 1998). We 
choose the NEW profile as the density profile of all halos: 



p{r) 



Ps 



(3) 



where the maximum circular velocity is achieved at 
Tmax — 2.16rs. Thc NEW profile is charaterized by the 
concentration parameter Cvir = Rvir/fs- Numerical simu- 
lations suggest several analytical functions with the mass 
of a halo and a redshift, Cvir = Cvir(M, z) (NEW; Bul- 
lock et al. 2001, hereafter BOl; Eke, Navarro & Stein- 
metz 2001; Maccio, Button & van den Bosch 2008). In 
our model, we use the model of BOl for the concen- 
tration parameter of a dark matter halo having mass 
Afvir at redshift z, which has collapsed at redshift Zc'. 
Cvir(M, z) = K{l + Zc)/{l + z), where BOl suggest K = 4 
since this agrees well with the result of their ACDM sim- 
ulations. This expression for Cvir(Af, z) indicates that 
massive halos are less concentrated than less massive ha- 
los and that a halo at low z is more concentrated than 

1/3 

that at high z. Since i?vir oc A^(^ (1 -I- z), an inner radius 
Vs is almost constant at all z. 

2.1.3. Orbital Evolution of Subhalos 

We track the orbital evolution of all accreted subha- 
los, taking into account the dynamical friction and tidal 
mass loss while orbiting in a host halo. Through these 
processes, subhalos eventually sink to the center of the 



host halo (termed "central merged" in ZB03), or lose the 
most of the mass ("tidal disrupted"), or keep the distinct 
structure from their accretion epoch to the present. We 
model these processes based on the ZB03 formalism. 

We define that a subhalo has a mass Mgat and outer ra- 
dius i?sat at thc accretion rcdshift of Zacc obtained by our 
merger tree, where i?sat is given by i?vir [equation 
The concentration parameter of a subhalo is calculated 
from thc median value as given by Afgat and Zacc at the 
accretion. Notice that i?sat and A/gat are varied due to 
mass loss while orbiting within a host halo. We update 
thc mass and concentration parameter of the host halo 
at each accretion event and continue to fix its density 
structure of the host halo while each orbit is integrated. 
In order to track thc orbit of each subhalo, wc assume 
the potential of the host halo to be spherically symmetric 
and static, so the orbit of a subhalo depends only on the 
radial distance from the center of the host halo. 

When it accretes onto the host halo, we assume that 
each subhalo has an initial orbital energy, which is con- 
sistent with the result of detailed N-body simulations by, 
e.g., Klypin et al. (1999). Namely, as described in ZB03, 
each subhalo has an initial orbital energy equal to the 
energy of a circular orbit at the radius i?circ = vRvin 
where i?vir is the virial radius at the accretion and rj is 
drawn randomly from a uniform distribution on thc inter- 
val [0.4, 0.75]. In this way, we define the absolute velocity 
of a subhalo at the accretion. Next we consider the ve- 
locity vector of a subhalo. We define the initial angular 
momentum of a subhalo at the accretion as J = eJcirc, 
where Jcirc is the angular momentum for a circular orbit 
of the same energy and e is a parameter, so-called "or- 
bital circularity" (LC93). Recent numerical simulations 
show that this orbital circularity is well fitted by beta 
distribution (Zentner et al. 2005): 
df{e) T{2a)^,_ 



-\i-er 



(4) 



de F2(a) 

where a = 2.22 provides a good fit to their ACDM simu- 
lation results with the mean (e) = 1/2 and the dispersion 
cr(e) = l/2(2a-hl)i/2. This form of the orbital circularity 
is also in agreement with the results obtained from cos- 
mological simulations by Benson (2005) and Khochfar & 
Burkert (2006). We define that the initial radial position 
of a subhalo at the accretion is i?circ and its orbit is given 
by dR/dt < in order to be initially infalling. For the 
calculation of its orbital decay by dynamical friction, we 
assume a subhalo as being represented by a point mass in 
thc gravitational potential of the host halo. Wc use the 
standard Chandrasekhar formula (Chandrasekhar 1943) 
for the evaluation of the frictional force, following the 
method shown in ZB03. 

The mass of a subhalo is stripped by the tidal force 
while orbiting in the potential of thc host halo. This 
process is approximately modeled as follows. First, the 
instantaneous tidal radius rt of the subhalo at each point 
along its orbit is calculated by the equation, in the limit 
that a subhalo is much smaller than a host halo, as (e.g., 
King 1962) 

3 Afgat(< rt)/A^host(< r) 3 

"^'-2 + c^V3/GA/irost(< r) - ainA^host(< r)/dlnr'' ' 

where r is the radial position of a subhalo in the host 
halo. A/host (< r) is the mass of the host halo within 
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r, Aisat(< Tt) is the mass of the subhalo within rt, to 
is the instantaneous angular speed of the orbiting sub- 
halo Lu = L/r'^, and L is the angular momentum of the 
subhalo. The process of the mass loss occuring out- 
side rt is calculated following Taylor & Babul (2001): 
we divide the orbit of the subhalo into discrete time 
steps of size 5t T, where T is the orbital time scale 
T = 27r/w, then at each time step, we remove the mass 
5m = Msat(> rt){St/T), where Msat(> '''*) is the mass 
of the subhalo outside rt. While a subhalo loses mass 
due to tidal force, its density profile is assumed to re- 
main unchanged within its outer radius i?sat, which is 
now defined as the radius within which the bound mass 
is contained in the form of an NFW profile. We continue 
to fix the scale radius of the subhalo, rf^*^, at the value 
calculated at the epoch of the accretion. 

Finally, we define whether a subhalo is tidally dis- 
rupted, centrally merged, or survived, following the cri- 
teria given in ZB03. We first set 7-^*^ — 2.16rf^ as the 
radius at which the initial circular-velocity profile of a 
subhalo reaches its maximum, and its mass A/gat (< ?'max) 
within rfjjax- We consider a subhalo to be "tidally dis- 
rupted" if its mass becomes less than Msat(< ?'max)- ^ 
subhalo is "centrally merged" with the host halo if its 
radial position, r, becomes smaller than rj^^ax- Also, 
the mass lost by tidal force modifies the circular-velocity 
profile of subhalos. To account for this effect, follow- 
ing ZB03, we determine whether or not the tidal ra- 
dius of each surviving subhalo has ever been less than 
r^ax- If the subhalo experienced so, we change the 
maximum circular velocity of the subhalo by V^^jax' ~ 
(M,*™'^VM™"'^V/3y™tiai^ ymitiai jg ^j^p maximum 

circular velocity of the initially defined subhalo, M^H^^ is 
its final mass, and M^a"*'^' is its initial mass at the ac- 
cretion. 

2.2. Model Checks 

Following ZB03, we have calculated the three examples 
for the orbits of subhalos. For all of these cases, we adopt 
the same initial condition of e = 0.5 and 77 = 0.5 for each 
subhalo's orbit, with the concentration and the mass of 
the host halo being Cvir = 6 and Mhost = 5 x 10^ ^Mq, 
respectively. The other parameters of the subhalos are 
given as follows: M° = lO^^M©, Cvir = 15; M°at = 
IO^Mq, Cvir = 7.5; A'/fat = 5 X IO^Mq, Cvh- = 15. The 
results are shown in Figure [Jl where the upper panel 
shows the evolution of radial positions of the subhalos 
and the lower panel shows the evolution of their masses 
of the subhalos. We set the start of calculation to be at 8 
Gyr ago, which is z ~ 1.14 for our cosmology. The plots 
in this figure arc in excellent agreement with Figure 3 in 
ZB03, indicating that our model successfully reproduces 
their results. 

We show, in Figure [21 the cumulative velocity function 
to demonstrate the statistical properties of our model. 
The plots correspond to an ensemble mean of 200 host 
halos with mass 1.4 x at z = 0, where these 

host halo systems have the maximum circular velocity 
Knax — 188 km s~^, the virial radius i?host — 285 kpc, 
and the concentration parameter Cvir — 13.4. We use 
^^min = IO^Mq. It is found that the cumulative velocity 
function of the survived subhalos until z = is in good 
agreement with Figure 4 of ZB03 and also the result of 
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Fig. 1. — Orbital evolution for three sets of the subhalo param- 
eters: Initially at t = 0, M°^^ = W^Mq, Cvir = 15 (solid Une); 
MO^j = 10**AfQ, Cvir = 7.5 (dashed line); M°^^ = 5 X IO^Mq, 
Cvir = 15 (dash-dotted line). The host halo is given by Cvir = 6 
and Mjiost = 5 X IO^^Mq. The top panel shows the radial evolu- 
tion in units of the initial radius as a function of time. The bottom 
panel shows the mass evolution in units of the initial mass as a 
function of time, where the vertical lines correspond to the end 
point when the subhalos are destroyed. 
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Fig. 2. — Cumulative velocity function of the survived and de- 
stroyed subhalos for the host mass 

Miost = 1.4 X lOi'*M0 at ^ = 0. 
These plots are an ensemble mean of 200 host halo systems. Each 
line is for the survived subhalos until the present epoch (thick solid 
line), tidally destroyed ones (dashed line), central- merged ones 
(thick dotted line), and all of the accreted subhalos as obtained 
from the merging history (dash-dotted line), respectively. 



the numerical simulation by Klypin et al. (1999). 

2.3. Method for Calculating Likely Disk Sizes 

A host halo contains numerous subhalos, among which 
there exist several very massive ones with mass of M ^ 
lO^°Af0. In this study, we assume that when such a 
massive subhalo having a comparable mass to a galactic 
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stellar disk of M ^ interacts with each other, 

then a stellar disk is made so thick due to dynamical 
heating that it no longer exists as a thin stellar disk like 
that of the Milky Way at the present epoch (e.g., Hayashi 
& Chiba 2006; Villalobos & Helmi 2008). More specifi- 
cally, we assume that if a subhalo having this large mass 
interacts with a galactic disk at any radial positions, then 
the latter is supposed to be destroyed by a subhalo. To 
account for this process in our model, we calculate, for 
the redshift interval oi z to z + dz, the minimum of peri- 
centric radii r^i^ for subhalos' orbital motions, where 
their mass is larger than a supposed mass of a galactic 
disk (Mth). We then suppose that this minimum ra- 
dius provides the upper limit for the radius of a galactic 
disk which can exist at z. It is thus postulated that any 
galactic disk having a size larger than this minimum of 
subhalos' pericentric radii will be destroyed by dynam- 
ical interaction with subhalos. In other words, the size 
of a galactic disk at z is limited by the orbits of massive 
subhalos in the redshift interval of 2; to z -I- c?z. 

We note that in the actual calculation of merger trees, 
there are only a few subhalos with mass M > Mth from 
z to z + dz, which makes it difficult to estimate the min- 
imum of pericentric radii in the interval z to z + in a 
statistically significant manner. Instead, we first consider 
the cumulative distribution from z to z = 0, for the mini- 
mum of pericentric radii for subhalos with M > Mth- We 
then convert this cumulative form with respect to z into 
the diiferential form from z to z + dz. Notice that in this 
calculation we exclude the subhalos which end up with 
"centrally merged" . Such subhalos sink to the center 
of the host halo, which results in major mergers with a 
galactic disk. Therefore, by excluding such subhalos, we 
can examine a dynamically quiet disk system that does 
not experience any major-merger events, thereby having 
a possibility to exist as a galactic disk. 

3. RESULTS 

3.1. Distribution of Pericentric Radii of Subhalos 

Here we show the fiducial case of a host halo with 
Mhost = 1.4 X 1O^^M0 at z = and the minimum mass of 
merger tree is Mmin = IO^Mq. Figure [3] -0 show the his- 
togram of the minimum of pericentric radii for orbiting 
subhalos with M > Mth from z to z = 0. In these fig- 
ures, we vary z from 0.3 to 1.8. FigurelHUH and[5]show, 
respectively, the cases of Mth = lO^^A/©, 6x 10^ Mq, and 
3 X IO^Mq, based on 200 merger trees. It is found that at 
fixed Mth, '"min at higher z is smaller than that at lower 
z. This can be explained as follows. At high z the virial 
radius of a host halo is small as suggested from equation 
([2]). Also the orbital time of subhalos which accreted at 
higher z is longer than those accreted at lower z. Thus 
the orbits of subhalos with longer orbital time decay more 
to inner radii since dynamical friction affects these sub- 
halos during a longer time, thereby yielding smaller r^i^. 
On the other hand, at fixed z, r^i^ for larger Mth is larger 
than that for smaller Mth- This may be because the ac- 
cretion of massive subhalos occur at preferentially low z, 
compared with less massive ones. 

Figure [5] shows the average of these histograms (plus, 
cross, and asterisk). The lines denote the fitting func- 
tions in terms of the least-squares method (without em- 
ploying a dispersion- weighted fitting), which are given 
as 22(1 + z)-0-54 for Mth = 10^°Mq, 17(1 + z)-"-^^ 



M„,=1.0x lO'^M,, 



M,h=1.0x10"'M„ 



0.3 
0.25 
0.2 

o 

g 0.15 

ro 

0.1 
0.05 


0.25 
0.2 

i 0.15 

ro 

^ 0.1 
0.05 


0.25 
c 0.2 
i 0.15 

S 0.1 

0.05 



ld3=ti:£LE£iL. 



I 



10 15 20 25 30 
■■min [kpcl 



10 15 20 25 30 
fmin [l<PCl 



Fig. 3. — The histogram of the minimum of pericentric radii for 
orbiting subhalos with the mass larger than Mjjj = 10^" Mq from 
the given redshift 2: to 2; = in each merger tree. In the upper 
right corner this 2 is depicted. We exclude the subhalos which end 
up with centrally merged. The grid of the histogram is 1 kpc. All 
the panels are based on the 200 merger trees. 
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Fig. 4.— The same as Figure [3] but for A/th = 6 X IO^A/q. 



for Mth = 6 X IO'^Mq, and 12(1 + z)-°-^^ for Mth = 
3 X 1O^M0. If a dispersion is used so that a fit- 
ting is weighted in favor of a low dispersion, we ob- 
tain somewhat shallower fitting curves: 19(1 -|- z)~°'^*, 
13(1 -h z)-°-27, and 11(1 -h z)-°-42, respectively. These 
plots show clearly that at fixed Mth, "rmin at high z is 
smaller than that at low z and at fixed z, rmin with mas- 
sive Mth is larger than that with less massive one. 

In addition to these cases with a threshold subhalo 
mass of Mth < lO^^M©, we additionally consider much 
more massive cases, i.e., Mth = 5 x lO^^Mo and lO^^Af© 
in view of recent studies (Stewart et al. 2008; Pur- 
cell et al. 2009) showing that such massive subhalos 
are most crucial in the survival of a thin disk. In Fig- 
ure [51 open and filled squares show the cases of these 
threshold masses, respectively. It is found that such 
massive subhalos, which do not end up with the cate- 
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Fig. 6. — Ensemble mean of the minimum of pericentric radii 
for orbiting subhalos with the mass larger than Mth = lO^^Mo 
(plus), 6 X IO^Mq (cross), and 3 X lO^M© (asterisk) from redshift 
2 to 2 = 0. Also, shown are the additional two cases with more 
massive subhalos, Afth = lO^^AfQ (filled square) and 5 X lO^^M© 
(open square). All the plots are based on the result of 200 merger 
trees. All the lines are the results of the least-squares fitting to 
these plots for the ensemble mean. 



gory of centrally merged ones, remain at large orbital 
distances of > 40 kpc and that the distribution of their 
pericentric distances stays nearly constant: the fitting 
functions are given as 44(1 + and 46(1 + 2;)+0 03 

for Mth = 5 X IO^^Mq and 10^^ Af©, respectively, for 
both with and without using a dispersion in the fitting 
procedure. Thus, a disk which is exempt from the central 
merging of such massive subhalos is likely able to have 
a large radius, unless the effects of less massive ones are 
significant. 

To derive the distribution from z to z + dz, we convert 
these fitting functions into differential forms as follows. 
We define N{z, i?sizc, Afdisk) as the distribution function 
of galactic disks with the size, i?sizc, and mass, Mdisk, at 
z. Considering the average of the size of galactic disks 




Fig. 7. — Redshift dependence of disk-size evolution in var- 
ious disk models and observed disks, where comparison is made 
by normalizing a disk size to be 1 at z = 0.1. Solid, long- 
dahsed, and short-dashed lines denote the fitting results (without 
considering the dispersion) to the average of the minimum peri- 
centric radius rmin in subhalos with mass which is larger than 
Afth = 10^° Mq, 6 X lO^ Mq, and 3 x 10^ Mq, respectively, from z 
to z + dz. Also shown are other model predictions: dotted line is 
the expected evolution of disk scale length from Mo et al. (1999) 
rs{z) /re{0.1) oc H~^/^{z), and dash-dotted line is the expected 
evolution of disk scale length from BS02 re{z) /re (0.1) oc 1 — 0.272, 
where is the effective radius. The points denote the evolution 
of an average disk scale length from observational result (Trujillo 
et al. 2006), re(2)/re(0.1). 



from the redshift z to z = 0, then we can write this as, 



/o°°/>(^'.«size,^disk)d^'di?^,e 



(6) 



Assuming that all disks are destroyed when subhalos with 
mass larger than the threshold mass Afth (comparable 
to a disk mass) interact with anywhere inside a galac- 
tic disk, (i?size) agrees with the minimum of pericentric 
radius for orbiting subhalos with mass larger than Afth 
from the redshift z to z = 0. We define the distribution 
of galactic disks with mass smaller than A/disk from z to 
z -f as i?sizc ■ We can write this as 



/o°° -^sizc^(^> -^sizc: MgiskldR'sv. 
/o°°^(~^^size,^disk)di?Le 



(7) 



With normalization N{z,R'^-^^^,Mdisk)dR'sizc — 1' 
obtain ^ 

(-Rsizc) = - / dz'R~{z'). (8) 
z Jo 



Assuming (i?sizc) 
tion. we obtain 



a(l -I- z)^ and employing differentia- 



i?size(^) = a(l + zf-\l + (1 + I3)z]. 



(9) 



In this way, we can convert the cumulative function over 
[0, z] into the differential function over [z, z + dz]. 

Figure [7] shows these differential functions normalized 
at z = 0.1, so that redshift dependence of size evolution is 
highlighted. It is evident that there is no significant dif- 
ference in redshift dependence between our models with 
different mass thresholds Mth- We also plot the predic- 
tion of the existing disk evolution models: Mo, Mao & 
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White (1999) model of r,{z)/rei0.1) cx H-^/^{z) (dot- 
ted line) and BS02 model of re{z)/re{Q.l) (x 1 - 0.27z 
(dash-dotted line), where re is the effective radius. Also 
shown (points with error bars) is the evolution of average 
disk scale length from observational results (Trujillo ct 
al. 2006). It follows that the redshift dependence of all 
the characteristic sizes is basically similar to each other, 
especially at z > 1. 

3.2. Comparison with Observed Disk Sizes 

Our interest in this work is to elucidate whether a 
galactic disk survives in the CDM model. For this pur- 
pose, we need to estimate the likely outer edge of an 
observed disk, which generally has an exponential light 
distribution. For a class of disk galaxies, this light dis- 
tribution does not extend to an arbitrarily large radius, 
but is radially truncated at the outer part. This radius is 
a so-called truncation radius, which can be regarded as 
the outer edge of the disk. Observations of nearby disk 
galaxies indicate that this truncation radius is related to 
the disk scale length as i?brcak ~ "^Rd — '^Rd (e.g. Kregel 
et al. 2002). In contrast, at high z few observations are 
available for the information of a truncation radius. We 
simply assume the relation -Rbrcak = 3.5i?(; at current 
epoch for such remote disk galaxies. 

We consider the Milky Way-type disk with mass 
lO^^M© and the scale length i?^ = 4 kpc for the present 
experiment, to obtain the evolution of the disk's outer 
edge. In Figure [H the points show the redshift evolution 
of a supposed disk size at fixed mass as derived from an 
observed disk scale length. Notice that we assume r^ is 
equal to Rd, since the effective radius is related to the 
scale length by a simple function. It follows that at low 
z the points are well below our theoretical upper limits, 
thereby suggesting that observed galactic disks at low z 
can survive in the CDM model. On the other hand, at 
high z there is a possibility that some of galactic disks 
may be destroyed by dynamical effects of subhalos, yield- 
ing some slight disagreement with observations. This 
possibility can be precluded by other effects occurred 
in a disk. For instance, when a galactic disk contains 
the large amount of cold interstellar gas, such gas which 
is unaffected by subhalos can reduce the thickness of a 
galactic disk. Also, at the truncation radius the mass 
density of a galactic disk are low enough that the effect 
of subhalos on a disk structure is minor. 

4. DISCUSSION AND CONCLUDING REMARKS 

We have investigated the dynamical effect of subhalos 
on the size evolution of a galactic disk, based on a semi- 
analytic method that includes approximated and empir- 
ical relations for the evolution of subhalos as obtained in 
detailed numerical simulations of the CDM model. This 
is motivated by our previous work (Hayashi & Chiba 
2006) that subhalos are able to induce the dynamical 
heating and tidal destruction of a galactic disk, thereby 
affecting the distribution of the size of a disk which does 
not suffer from the interaction with subhalos; such a disk 
ought to be smaller than the minimum of pericentric or- 
bital radii of subhalos. We have found that the upper 
limit for the size of a galactic disk with large mass is 
larger than that with small mass. Also, with fixed mass, 
the upper limit for the size of a galactic disk at low z is 
larger than that at high z. 
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Fig. 8. — The same as Figure [7] but without being normalized. 
We assume a Milky Way type disk, where re (0.1) = 4 kpc and the 
size of a disk is 3.5re. 



Stellar-population analyses of present-day galaxies 
have revealed the so-called 'downsizing' evolution, where 
the stars in more massive galaxies tend to have formed 
earlier and over a shorter time-span (Thomas et al. 2005; 
MacArthur et al. 2004). On the other hand, in the CDM 
model halos are formed hierarchically in a bottom-up 
manner, contrary to the downsizing evolution. In view 
of this general downsizing evolution of galaxies, our result 
that a galactic disk with larger size and mass is available 
only at lower z may conflict with observations. However, 
it is noted that this evolution of a disk is limited by the 
pericentric orbital radius of subhalos, i.e., a disk which 
has never been affected by subhalos, whereas a galactic 
disk having initially a large radius can indeed be dy- 
namically heated and/or tidally dcstructed by subhalos. 
Even after such a destruction incident of a thin stellar 
disk, residual gas or infalling fresh gas allows to form a 
new disk component. Indeed according to recent numer- 
ical simulations, infalling subhalos do not fully destroy a 
galactic disk, where some fraction of the thin disk com- 
ponent are survived (Kazantzidis et al. 2008; Villalobos 
& Helmi 2008). We thus regard that old stellar com- 
ponents in massive disk galaxies may be the relic of a 
first-existing galactic disk. In particular, recent numeri- 
cal simulations have demonstrated that a thick disk can 
be formed by the dynamical interaction between such a 
first-existing disk and subhalos (Hayashi & Chiba 2006; 
Kazantzidis et al. 2008; Villalobos & Helmi 2008). 

Several models have been proposed for the evolution 
of a disk scale length (e.g., MMW98; BS02). However, 
such models do not explicitly take into account the fi- 
nite extent of a galactic disk. For instance, the MMW98 
model assumes that stars form from gas instantaneously, 
while maintaining an exponential density profile of a disk. 
Since an exponential disk does not have an outer limit in 
its profile, this model cannot provide the size evolution 
of a galactic disk. On the other hand, the BS02 model 
adopts the present state of a disk galaxy as an initial 
condition for its backward approach in time. This model 
also assumes an exponential stellar disk at present, so it 
is unclear how the disk size evolves; we need to take into 
account explicitly the size evolution of a disk in its model 
for the formation of an exponential density profile. 
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As such a model, we consider a so-called viscous disk 
model (Lin & Pringle 1987; Yoshii & Sommer-Larsen 
1989). This model shows that the exponential distri- 
bution of stars in a galactic disk is a natural product 
of angular momentum redistribution caused by viscosity. 
Yoshii & Sommer-Larsen (1989) found that if the con- 
dition that star formation time scale being comparable 
to viscous time scale is maintained, then the final stel- 
lar distribution shows an exponential profile irrespective 
of the specific form of the initial gas profile. However, 
since stars are assumed to form continuously as long as 
gas is still remained, where fresh gas is always replen- 
ished by infalling gas, it is evident in this model that no 
outer limit comes out in a galactic disk. We thus need 
to consider, for instance, the star formation threshold 
in this model, to obtain the size evolution of a galactic 
disk, where stars are formed when gas density is beyond 
some threshold value. Such a threshold in gas density 
may be defined by a so-called Toomre criterion for gravi- 
tational instability of an axisymmetric disk and/or some 
external feedback such as a UV background radiation, 
although an explicit account of these processes is yet to 
be explored in a viscous disk model. 

It is worth noting that in our models presented here, 
a couple of other effects, e.g., gas-dynamical response of 
a disk and tidal force of a disk on subhalos are ignored. 
When we consider cold gas components in a disk, the 
stronger binding energy of such a disk than that without 
gas may be less affected by subhalos. Also, the tidal force 
of a disk on subhalos may make the mass loss of subhalos 
more significant. Both of these effects help the survival 
of a galactic disk. If we assume that the orbital circu- 
larity of subhalos is more tangential than that adopted 
here, the number of subhalos having nearly circular or- 
bits becomes larger. Then subhalos can be more easily 
survived at larger orbital radii, so the size of a galac- 
tic disk limited by such subhalos may be larger. This 
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therefore suggests that it is essential to understand the 
realistic dynamical properties of subhalos, based on the 
CDM model, in order to set tighter constraints on the 
evolution of a disk galaxy like the Milky Way. 

We conclude that subhalos tightly regulate the struc- 
ture and evolution of a galactic disk. The effects of mas- 
sive subhalos on a galactic disk are significant compared 
to less massive ones, so it is important to consider the dy- 
namics and properties of such subhalos. Since at high z 
subhalos can sink to the central region of a host halo, the 
dynamical effects of subhalos on a galactic disk is more 
efficient at high z than at low z. Then, a galactic disk 
at present might experience the interaction with subha- 
los one or more times during its evolutionary process, 
especially at high z. We postulate that a thick disk is a 
relic of a pre-existing galactic disk after these interaction 
events. 

It is also important to construct the new formation 
model of a galactic disk taking into account its size evo- 
lution and dynamical effects of subhalos, so that one can 
directly compare the model results with observed disks 
at high z. In this respect, it is useful to carry out high- 
resolution hydrodynamical simulations and predict de- 
tailed properties of galactic disks at high z. Further- 
more, more detailed observations of extragalactic thick 
disks, in particular their kinematics, are needed to find 
the evidence for the interaction between subhalos and a 
galactic disk (e.g., Herrmann, CiarduUo, & Sigurdsson 
2009). 
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